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ABSTRACT 

m , 

, The spectrometer SPI, one of the two main instruments of the INTEGRAL spacecraft, has strong capabihties 
■ in the field of Gamma-Ray Burst (GRB) detections. In its 16° field of view (FoV) SPI is able to trigger and 
I to localize GRBs. With its large anticoincidence shield (ACS) of 512 kg of BGO crystals SPI is able to detect 
• GRBs quasi omnidirectionally with a very high sensitivity . The ACS GRB alerts will provide GRB arrival 
D ' times with high accuracy but with no or very rough positional information. The expected GRB detection rate in 
SPI's FoV will be one per month and for the ACS around 300 per year. At MPE two SPI software contributions 
I to the real-time INTEGRAL burst-alert system (IBAS) at the INTEGRAL science data centre ISDC have been 
developed. The SPI- ACS branch of IBAS will produce burst alerts and light-curves with 50 ms resolution. It is 
planned to use ACS burst alerts in the 3'''^ interplanetary network. The SPI-FoV branch of IBAS is currently 
under development at MPE. The system is using the energy and timing information of single and multiple 
^ir) I events detected by the Germanium-camera of SPI. Using the imaging algorithm developed at the University of 
Birmingham the system is expected to locate strong bursts with an accuracy of better than 1°. 
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1. INTRODUCTION 



Q^' The launch of INTEGRAL, ESA's next mission devoted to 7-ray research, is expected soon.-'^'^ The scheduled 

' I launch date is October 17, 2002. A Russian PROTON rocket will carry the spacecraft from the Cosmodrome 

^ ■ in Baikonur, Kazakhstan on to its higly eccentric orbit with a revolution period around the Earth of 3 days. It 

I has a perigee height of lO'OOO km and an apogee height of 152'600 km with an inclination of 51.6". This orbit 

■ allows long unbroken observations outside the radiation belts. 

After a period of about 3 months, which is dedicated to the commissioning phase and the following inorbit 
k>( \ calibration and the scientific performance validation , INTEGRAL will start its scientific observation program.^' ^ 
5_j ■ The detection and investigation of cosmic gamma-ray bursts (GRBs) is one of the important scientific topics of 
C3 : the INTEGRAL mission. 

GRBs are still one of the most fascinating research topics in astrophysics. Until now, since their discovery 35 
years ago by the Vela satellites in 1967,^ this phenomenon is still not totally understood and explained. The first 
major breakthrough in this field was obtained with the BATSE detectors^ on NASA's Gamma Ray Observatory 
GRO mission. In the ten years of the GRO mission around 3000 bursts were registered, which showed an 
isotropic distribution over the entire sky, but with a deficiency of weak bursts. It was the Italian/Dutch satellite 
BeppoSAX^ which revealed the cosmological nature of GRBs with the identification of the first X-ray afterglow 
in 1997,^ which triggered the first successful follow-up observation at optical wavelengths.* This finally ruled 
out the Galactic population models. In the meantime an optical afterglow was measured for 25 GRBs, all of 
them showing redshifts between z=0.36 and z=4.5.^ Many new missions, especially built for burst searches or 
equipped with instruments for the investigation of GRBs are planned, will be launched soon^^ or are already 
in orbit. 
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^Prom all instruments of the INTEGRAL payload, the two main instruments, the spectrometer SPI^^~^^ 
and the imager IBIS,"'^^ and the two monitors, JEM-X-'^'^ and OMC,""^^ contributions to the science of GRBs are 
expected. Especially the sensitivity of INTEGRAL'S imager and spectrometer above 20 keV will distinguish 
them from other missions. Although INTEGRAL is not equipped with an on-board burst alert system it has 
strong capabilites to detect bursts thanks to the on-gromid-bascd INTEGRAL burst-alert system IBAS, which 
was proposed in 1997^^ after the success of the BeppoSAX satellite in rapidly localizing GRBs. 

IBAS is an automatic software system for the near-real time distribution of GRB alerts detected by IN- 
TEGRAL to the scientific community. The telemetry stream of the satellite, which is linked via the mission 
operation center (MOC) at Darmstadt, Germany, to the INTEGRAL science data center (ISDC)^° at Versoix, 
Switzerland, will be monitored continuously for the occurrence of GRBs. The expected transmission time of 
the telemetry from the satellite to ISDC will be approximately 30 s. Additionally the processing time at ISDC 
and the time needed to distribute the data to the clients has to be added. In total an observer could expect 
the first alert message within less than 1 min after the onset of the burst. Thus, in the case of a long-duration 
burst, follow-up observations are possible even during the active phase of the burst. 

The principal structure of the IBAS system is displayed in figure 2, showing the programs dedicated to the 
burst search for the two main instruments IBIS and SPI (two each), for the anticoincidence system ACS of SPI 
and for JEMX. This report will focus on the three software packages written for the SPI and SPI/ACS. An 
overview of IBAS with details on IBIS are shown in several proceeding reports. ^^"^^ 

2. THE SPECTROMETER SPI 

The aim of the spectrometer SPI is to perform high-resolution spectroscopy in the energy range between 20 
keV and 8 MeV.^^~^^ The imaging capability will be good, but it is exceeded by that of the imager IBIS 
which complements SPI with its high imaging resolution, but with less spectroscopic resolving power. So the 
expected yield and discoveries in the field of GRB research will be of a different kind for each of the four 
INTEGRAL instruments. With SPI new discoveries in the spectroscopic regime are possible, supposing that 
previously- unkown spectral features exist. IBIS will provide more precise GRB locations to the community, 
which is important for the observation of GRB afterglows. The minor sensitivity of SPI for GRB detection is 
compensated by its larger FoV. So the expected rate/year of detected GRBs is about the same for IBIS and 
SPI. All these quantities arc summarized in table 1. Also the two monitor instruments will help to augment the 
scientific output for GRBs of the INTEGRAL mission. 

Table 1. This table shows the important performance parameters for the burst detection by the two main INTEGRAL 
instruments. PSLA stands for "point-source location accuarcy". 
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Fig. 1 shows a drawing of the spectrometer SPI. The camera of SPI, which consists of 19 cooled high-purity 
germanium detectors, residing in a cryostat. is shielded on the side walls and rear side by a large anticoincidence 
shield (ACS). The field of view (FoV) of the camera is defined by the upper opening of the ACS. A detailed 
description of the ACS can be found in section 4. The imaging capability of the instrument is attained by a 




Figure 1. INTEGRAL Spectrometer SPI. 

passive-coded mask mounted 1.7 m above the camera. Below the mask a plastic scintillator anticoincidence 
(PSAC) allows a reduction of the 511 keV background, which is mainly generated by particle interactions in the 
mask. By observing a series of nearby pointing directions around the source(" dithering" ) the imaging capability 
of SPI is improved by reducing ambiguity effects. For the normal mode of operation two dithering patterns 
are planned, a quadratic 5x5 and a hexagonal with 7 pointings. The pointing on one grid point will last for 
about 20-30 min with a short 5 min slew between succeeding pointings. Due to the short duration of the bursts, 
dithering will not improve the imaging of SPI for bursts. But during its short duration a burst will be in most 
cases the brightest 7-ray source for SPI in the whole FoV, so imaging with SPI is still possible. For that it is 
required that the true source position corresponds to the point in the image with the most-significant detection. 
The possibility of observing a burst occurring during a slew, has still to be investigated. Despite the modest 
angular resolution of SPI, which is in the order of 2° , it is possible to locate the direction of bursts down to a 
few arcminutes (see chapter 8.1). 

3. BURST DETECTION IN SPI'S FIELD OF VIEW 

The triggering algorithm which will be used for the burst search in SPI's FoV will be similiar to that used 
inside IBAS for the IBIS-ISGRI instrument. The first algorithm is based only on imaging. It looks for the 
appearance of new sources by comparing the actual image of the sky with those obtained previously. This 



algorithm is more sensitive to slowly rising bursts, since the imaging routine needs some time. The second is 
expected to be faster, because it monitors the overall event rate in all detectors together as function of time. In 
the case of a trigger the second algorithm uses the imaging algorithm for verification. It will test if the count- 
rate increase is connected with the appearance of a point-like source in the FoV, thus avoiding false triggers 
from background variations and other instrumental eflfects. Both algorithms can run in parallel with multiple 
instances, using different sets of parameters like energy interval or integration time. 

Every interaction of a 7-ray in one of the Ge-detectors is down linked event-by-event, wrapped into packets 
within the telemetry stream. The event-by-event packets contain information on the detector ID which was hit, 
the energy channel of the recorded pulse-height signal and the time of interaction, with an accuracy of about 
100 IIS. Two kinds of interactions of 7-rays in the Ge-camera have to be distinguished, in the first case a 7-ray 
only interacts within one of the Ge-detectors (single detector events), in the second case the 7-ray is interacting 
with more than one detector (multiple detector events). The condition for the occurrence of a multiple-detector 
event is defined by a coincidence window inside the digital front-end electronics (DFEE) of SPI. Both types of 
multiplicity will be downlinked with the science telemetry packets, and are used both for the GRB detection. 

Additional information on the detector interaction is available through the science house keeping data (HK- 

data). The rate for each detector is recorded on a 1 s basis inside the DFEE and sent to ground, but without 
energy information. The advantage of these data is that they are always available, even in the case of telemetry 
limitations. For the burst case a 3 Mbit buffer is used on board, and the read out of the buffer could last 
several minutes, depending on SPI's actual polling-sequence table (PST). Another application of the science 
HK Ge-rates could be the case that SPI is working in emergency mode. In this mode only the spectra and the 
multiple-detector events are included in the telemetry. 

4. THE ANTICOINCIDENCE SHIELD OF SPI 

One important part of SPI is its large anticoincidence shield ACS which provides a large effective area for the 
detection of bursts, but unfortunately with no or small positional information. But adding ACS to the 3'''^ 
IPN will avoid this weakness. 

The parts of the ACS are also shown in Fig. 1. The ACS consists of 91 BGO crystals which are arranged in 4 
subunits (Fig. 1). The units of the upper veto shield (UVS) are composed of the upper-coUimator ring (UCR), 
the lower-collimator ring (LCR) and the side-shield assembly (SSA), each containing 18 crystals which are 
arranged hexagonally around the cylindrical axis of SPI. The lower veto shield (LVS), consisting of 36 crystals, 
is assembled as a hexagonal shell. The thickness of the crystals increases from 16 mm at the top (UCR) to 50 
mm at the bottom (LVS). The total mass of BGO used for the ACS is 512 kg resulting in the obvious use of 
the ACS as a burst monitor. 

Each BGO crystal of the ACS (with one exception) is viewed by two photomultipliers (PMTs). Due to the 
redundancy concept used for the ACS, each of the 91 front-end electronic boxes (FEEs) sums the anode signals 
of two PMTs, which view different BGO crystals (in most cases neighbouring crystals). This cross strapping of 
FEEs and BGO units avoids the loss of a complete BGO crystal if a single PMT or FEE fails. It emerges that 
a disadvantage of this method is an uncertainty in the energy-threshold value of individual FEEs, caused by 
different light yields of neighbouring BGO-c;rystals and different PMT properties like quantum efficiency and 
amplification. A result of this is that the threshold extends over a wide energy range and is not at all sharp. 
The actual energy-threshold settings of the ACS will be the result of a tradeoff between background reduction 
and deadtime for the SPI camera. 

5. THE ACS AS GRB MONITOR 

The main task of the ACS as a detector is the generation of a veto signal for the camera to suppress charged 
particles and 7-rays coming from outside the FoV. But the signals from the ACS can also be used for scientific 
purposes. The ACS housekeeping (HK) data include the values of the overall veto counter of the veto control 
unit (VCU) and the individual ratemeter values of each FEE. Both kinds of HK data are suitable for burst 
detection. The count rate of the overall veto counter (ORed veto signals of all 91 FEEs) is sampled every 50 



ms. A packet, containing 160 consecutive count rates, will be transmitted every 8 s to ground. If no gap in 
the telemetry stream occurs one could have a continuous ACS veto-rate hght curve with 50 ms binning. The 
measurement time of the individual FEE ratemeter can be adjusted between 0.1 and 2 s. All 91 FEEs are 
read out successively in groups of 8 FEEs every 8 s. The readout of all 91 ratemeter values thus needs 96 s. 
In contrast to the VCU overall veto counter the individual ratemeter values do not yield a continuous stream. 
Additionally the values of different FEE groups are shifted by a time interval of 8 s. So it is very difficult to 
derive the burst-arrival direction from these individual counting rates. 

6. ACS AS PART OF THE IPN 

The 3'''^ interplanetary networks (IPN) of 7-ray detectors is in operation since the launch of the Ulysses spacecraft 
in late 1990.^^ Many spacecrafts and experiments have been included since that time. During the first year 
of the INTEGRAL mission the IPN will consist of ULYSSES, Mars Odyssey 2001, Konus-WIND, HETE-II, 
RHESSI and of course INTEGRAL (SPl/ACS).^'^ The network will have an excellent configuration during this 
time, due to the large spacecraft separations between Earth, Ulysses and Mars. 

An IPN localization of a GRB is obtained by comparing the burst's arrivel times of two spacecrafts {ti — ^2 = 
ST), which are separated by a long distance D. If the exact positions of the spacecrafts are known, one can 
localize the burst arrival direction to an annulus which spans from the spaceraft interconnection line with the 
half-angle 8 given by the following expression: cos0 = c • 5T/D (c: speed of light). The important quantity, 
the annulus width is given by AO = c • a (ST) / D sin O , which is mainly influenced by the uncertainty of the 
arrival-time-difference cr{ST) and the spacecraft separation D. In the above-mentioned configuration of the IPN, 
with INTEGRAL as near-Earth counterpart, annuUi widths in the order of about 1' could be obtained. This is 
due to the high time resolution of the ACS of 50 ms and the spacecraft separations of about 1500 lightseconds. 

With a third spacecraft a second annuUus can be determined, which has two intersections with the first one. 
This ambiguity can only be solved if one of the spacecrafts provides additionally a rough positional information. 
With four or more spacecrafts the burst location is known unambiguously. 

7. THE SPI AND ACS IBAS SOFTWARE MODULES 

Figure 2 gives a simplified overview on the structure of the whole IBAS system. After the near-real-time data 
receipt of the telemetry from MOC, there are several GRB detection programs running in parallel and searching 
in the telemetry stream for the occurrence of a burst. Three of them are responsible for the spectrometer SPI, 
there are the two programs monitoring the camera events (ibas_smonitor_inig and ibas_sinonitor_rate) and 
the monitor program for the ACS overall veto rate (ibas_acsinonitor). The other half of the branches monitor 
the science data of IBIS ISGRI (ibas_imonitor_img and ibas_imonitor_rate) and JEMX (ibas_j monitor). 
The ibasalertd program receives the GRB triggers from the individual branches and performs the final alert 
verification and distribution of the alert messages to subscribed users. 

All IBAS processes are multithreaded applications and run as daemon processes. So IBAS is able to perform 

several subtasks at the same time, with the advantage for the burst search to monitor the telemetry in different 
energy bands and with different time binning. Each of the GRB detection processes uses the frontend library 
to read out the telemetry packets provided by the data receipt. Then the required instrument data (Ge 
single/multiple events, Ge detector count rates, ACS veto rates) arc decoded from the telemetry packets and 
put into a time-sorted queue, which is large enough to hold several minutes of data. Several trigger-algorithm 
threads can access the data in the queues, in order to detect bursts with different sets of parameter values. 

The next two subsections present the ACS and SPI-FoV GRB detection programs in more detail. We 
begin with the ACS burst-alert system, because this software package has already been delivered to ISDC. The 
SPl-FoV software is still under development and the delivery of the first version is expected soon. 
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Figure 2. Structure of the IBAS. A more detailed description can be found in the "INTEGRAL Burst Alert System 
Architectural Design Document" of the ISDC (IBAS- ADD). 

7.1. The ACS Burst- Alert System 

The SPI/ACS Burst Alert System ibas_acsmonitor (an earlier version of the program was named sacs_grbtd, 

"SACS-BAS"^''') is one branch of IBAS (sec Fig. 2). The trigger algorithm used for ibas_acsmonitor looks 
for a significant excess with respect to a running average, comparable to the trigger algorithm used for other 
spacecrafts (e.g. ULYSSES). Up to now ibas_acsmonitor reads only the overall- veto counter values. But it is 
planned also to include the read out of ratemeter values of individual FEEs. This will allow a rough estimation 
of the GRB arrival direction. An accuracy of about 10° - 20° would be enough to distinguish between the two 
arrival-cone intersections of the interplanetary network (IPN). The ibas_acsinonitor trigger algorithm is being 
tested with generated telemetry data of simulated burst data plus background. After launch all parameters for 
the trigger algorithm and verify criteria will be optimised. 

Fig. 3 shows the thread structure of the ibas_acsmonitor program. The GRB detection threads (" detector" - 
threads) search, with different parameter settings (i.e. timescales, in order to be able to trigger on bursts with 
different temporal behaviour), for GRBs in SPI ACS light curves. If a significant increase is detected, an 
analysis (trigger verification) is performed by the temporarily created verify threads, to gain for the supression 
of fault triggers. In the case of a positive verification result, an alert is sent to the alert management system 
(ibasalertd). The number of GRB detection threads is likewise configurable via a parameter. Beside the 
GRB detection and verify threads, ibas_acsmonitors main program thread creates several other threads to 
perform different tasks: The signal handler thread processes all signals delivered to the ibas_acsinonitor 
process. The data receipt interface thread manages the telemetry data transfer, connects (and reconnects, in 
case of connection breakdown) to the near-real-time telemetry daemon rttmd process via TCP/IP socket. It 
reads the telemetry data packets, processes them and puts data into the local queue. The main detector thread 
initializes GRB detection threads, including the read in of the program parameters. Some checks are performed 
to guarantee correct variable ranges. The main verify thread watches the results of the verify threads and makes 
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Figure 3. Thread structure of the ibas_acsmonitor. 



the final decision about sending a burst alert to ibasalertd. If a burst is recognized and an alert must be sent, 
this thread sets all alert variables defined in alert format. Fig. 4 gives a rough overview of ibas_acsmonitors 
program sequence. 

Up to now several verification criteria can be used by the verify threads which are all ORed at the end. The 
idea of the verification tests/criteria is to use characteristic features of the burst light curve to distinguish them 
from other disturbances. But the practicability of each criterium has to be tested after launch. All the criteria 
are configurable via the parameter file. Due to the lack of knowledge of the the burst direction, the distinction 
between false and true alert is difficult (e.g. in the case of solar flares). 

The IBAS alert daemon ibasalertd has the task to distribute the GRB alert messages to a group of 
recipients (external users) which arc subscribed to the IBAS-alert distribution list. IBAS alerts arc delivered as 
alert packets via Internet using UDP transport protocol, which means alerts can be delivered to the recipients 
within 1 s. Different types of alert packets are available, with one of them especially adapted for the SPI/ACS 
alerts, containg information on the trigger time in universal time (UT), the significance of the GRB detection, 
the spacecraft position and its attitude (these last two informations are currently not containd in the SPIACS 
alert packets. In the next release of the ibasalertd this information should be included!). 

Especially for the IPN the burst time history (~ 100 s) is important for the alignment of the light curves 

obtained from different spacecrafts. For this purpose the time history together with the pre-trigger time-history 
(~ 5 s) will be transmitted to the IPN (up to now it is not decided how: via alert package or by making them 
available on the WWW after a short time). It is important for the IPN to know the burst-arrival time with 
millisecond accuracy. As already shown in^^ this is possible for the ACS overall-counter values. 

7.2. The SPI FoV Burst-Alert System 

The programs for the detection of GRBs inside SPFs FoV comprise the two modules ibas_smonitor_img and 
ibas_smonitor jrate (see Fig. 2). Both are currently developed at our institute, and the main structure of the 
programs has been taken over from the ACS program ibas_acsinonitor with a subsequent adaption to the new 
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Figure 4. Main steps performed in ibas_acsmonitor program. 

tasks and different data structure. The source-detection imaging algorithm/code was delivered by P. Connell, 
University of Birmingham. The code of the imaging module is part of the ibas_smonitor_ing ( jrate) packet 
and the imaging routine is called by the GRB detection threads. 

^From the telemetry packets the information on the 7-ray interaction time, the detector ID, and energy 
channel are extracted (see also chapter 3). Up to now single and multiple events are analysed. This is sufficient 
to do burst imaging in the instrument system. If a burst is detected and its position in the instrument system 
is determined, the sky coordinates can be computed subsequently from the pointing information. The energy 
channel is used to perform a basic energy selection, but an accurate channel-to-energy calibration is not required. 
The response will be prc-computcd for the standard energy channel selections. The data will be accumulated 
for a time determined by the parameter table and the burst search performed. This is repeated at least at the 8 
s rate of arrival of new data. However the accumulation time can be longer or shorter than this interval; several 
accumulation times can run in parallel to optimize detection of long and short bursts. The parameters should 
be as similar as possible to those for IBIS. They include: accumulation time, significance level, and energy 
channel range. 

For the delivery of the SPI-FoV alerts several burst packets can be distributed by ibasalertd. Directly 
after the GRB detection a packet containing a first rough estimate of the burst position is sent. Later packets 
with a refined position will be distributed. 

8. SENSITIVITY ESTIMATION 
8.1. GRBs detected in SPI's FoV 

Simulations and calculations have been performed^^ in order to investigate SPI's capabilities for the burst 
detection in its field of view. A GEANT Monte-Carlo software, adapted for SPI and a model of the instrument 
including a realistic distribution of the matter in the mask, veto-shield and detector assembly was used. By using 
a strong burst from the BATSE catalogue as input, other bursts were generated by scaling flux and duration. 
For a range of GRB arrival directions the detector count rates were simulated and a random background added. 




Figure 5. Effective area of the ACS in cm^ as a function of GRB arrival angle and energy. The SPI viewing direction 
is 0°. The IBIS instrument occupies the space to the left, leading to a reduced effective area. 

Finally these rates were analyzed with the independent software system CAPTIF^*^ with the result that SPI 
exhibits a relatively good performance for the detection of GRBs. For bursts with a fluence of 1.1 x 10^^ crg/cm^ 
and a duration of 1 s (only accounting for single detector events in the energy range from 20 to 300 keV) the 
positioning errors inside the FoV (< 9°) will be between 1.5' and 20' and for bursts with a larger offset (< 15°) 
the position error will still be less than 1°. Problems arising from ambiguities and systematic errors in the 
determination of the position are more serious towards the edge of the FoV. In summary GRB source locations 
to 10' are possible over a wide field of view. For weak bursts the detection could be problematic, the imperfect 
imaging together with random fluctuations can lead to a spurious peak appearing to be stronger than the true 
one. The best-estimate position can then be far from the correct location. The ambiguity problem can occur 
for bursts with a strength of less than 10~^ erg cm~^. 

As described in section 3, one method for detecting bursts in the SPI data is the monitoring of the overall 
event rate in all detectors together as function of time. The minimum-size burst detectable for a backgound 
rate of 0.25 cts cm~^ s~^ in the 50 to 300 keV regime for a 1 s burst, is about 0.2 photons, corresponding 
to 4 X 10~^ erg cm~^. This is about the same as the weakest bursts detected by BATSE. 

The rate at which GRBs are detected by an instrument depends on the FoV as well as on the sensitivity. 

SPI can detect bursts at off-axis angles up to at least 15°, so the useful FoV is quite large, 0.2 sr or 680 square 
degrees and consequently the expected detection rate of about 1 burst per month is comparable for IBIS and 
SPI. 

8.2. ACS effective area and expected GRB rate 

An estimation of the expected rate of GRBs detected by ACS has already been given.^^ For an effective area of 
about 3000 cm^, an ACS background rate between ~ 80000 cts/s and ~ 160000 cts/s, an ACS threshold of 80 
keV and a time binning of 50 ms the minimum detectable energy flux is 2 — 2.8 x 10~^ . For a time binning 
of 1 s the sensitivity improves to 5 x 10~ Tjfjfrj- Using the logN-logP distribution, measured by BATSE and 



PVO^^ one can derive the number of bursts which will be observed with the ACS in one year. The resulting 
values are ~ 50 bursts/year for 50 ms integration time and 280 bursts/year for 1 s integration time. The 
exact response of the ACS depends on the burst arrival angle, due both to projection effects and to shielding by 
neighboring instruments and by the spacecraft structure (sec Fig. 5). These effects together with occultation 
by the Earth on the ACS-"FoV" was not taken into account by the above-mentioned estimation of Lichti et 
al. 2000.^* The burst intensity can be determined once the infalling direction is known. This is possible only 
via Monte-Carlo simulations using the INTEGRAL mass model. Similar simulations have been performed by 
P. Jean et al.^^ in order to determine the sensitivity of the ACS for detection of novae. 
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